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I. INTRODUCTION 



Nonzero ground state disorder and associated entropy, 5*0 7^ 0, is an important subject in statistical 
mechanics. One physical example is provided by ice, for which the residual molar entropy is Sq — 0.82 ± 
0.05 cal/(K-mole), i.e., Sq/R — 0.41 ± 0.03, where R — NAvogkn Indeed, residual entropy at low 

temperatures has been observed in a number of substances, including nitrous oxide, NO, carbon monoxide, 
CO, and FCIO3 (a comprehensive review is given in Ref. [Q). In these examples, the entropy occurs without 
frustration. In magnetic systems, two examples are provided by the Ising antiferromagnet on the triangular 
and kagome lattices here, the ground state entropy does involve frustration. One should also mention 
the recent experimental observation of compounds whose behavior can be modelled by quantum Heisenberg 
antiferromagnets on a kagome lattice, including SrCr8_2;Ga4+2;Oi9, (SCGO(x)) and deuteronium jarosite, 
(D30)Fe3(S04)2(OD)6 [||. The quantmn Heisenberg antiferromagnet has a disordered ground state with 
finite entropy and frustration on the kagome lattice pO[ |. A particularly simple model exhibiting ground 
state entropy without the complication of frustration is the q-state Potts antiferromagnet (AF) 
lattice A, for q > x(A)j where xi-^) denotes the minimum number of colors necessary to color the vertices of 
the lattice such that no two adjacent vertices have the same color. As is already evident from the foregoing, 
this model also has a deep connection with graph theory since the zero-temperature partition function of the 
above-mentioned g-state Potts antiferromagnet on a lattice A satisfies Z{A,q,T — 0)paf — P{-A,q), where 
P{G, q) is the chromatic polynomial expressing the number of ways of coloring the vertices of a graph G ^ 
with q colors such that no two adjacent vertices (connected by a bond of the graph) have the same color 
|p3|-|l8|. Hence, the ground state entropy per site is given by S'o/fcs = In W^(A, q), where VF(A, q), the ground 
state degeneracy per site, is 

W{A,q) - lim P(A„,g)l/" (1.1) 

n — *oo 

Here, A„ denotes an n-vertex lattice of type A (square, triangular, etc.), with appropriate (e.g., free) boundary 
conditions. 

Since nontrivial exact solutions for this function are known in only a very few cases (square lattice for 
q — 3 triangular lattice ||2^], and kagome lattice for q ~ 3 ||l^), it is important to exploit and extend 
general approximate methods that can be applied to all cases. Such methods include rigorous upper and 
lower bounds, large-g series expansions, and Monte Carlo measurements. Recently, we studied the ground 
state entropy in antiferromagnetic Potts models on various lattices and obtained further results with these 



^We recall the strict mathematical definition of a graph as a collection of vertices and bonds connecting (all or a 
subset of the) vertices, such that there are (i) no loops, i.e., bonds connecting a vertex to itself, (ii) no multiple bonds 
connecting a given pair of vertices, and (iii) the total number of vertices is finite. Thus, an infinite lattice is not a 
graph, but a limit thereof. A graph is termed empty if it has no bonds connecting the vertices. 
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three methods |2l|-|24[|; see also Ref. [Q. As P{G,q) is a polynomial, q may be continued from integer to 
complex values, which we do here. 

One can gain valuable insight into the values of the ground state degeneracy of the Potts antifcrromagnet 
W{A,q) for different 2D lattices A and values of q by studying this function on infinitely long strips of 
increasing widths. Accordingly, in this paper, we calculate P{{Gs)rm q) for a family of strip graphs {Gs)m 
which are comprised of m repetitions of a subgraph unit H attached to an initial subgraph / on one end 
(which, by convention, we take to be the right end). Symbolically, we denote such a strip graph as 

m 

{GsU^{[[H)I (1.2) 

£=1 

for m > 1, with {Gs)q = The initial subgraph / may itself be equal to H. To carry out these calculations, 
we construct a generating function T{Gs, x) that yields the chromatic polynomial P{{Gs)rm q) as the 
coefficient of the term a;™ in its Taylor series expansion about a; = 0, where x is a formal auxiliary variable. 
We then use this generating function to compute the asymptotic limiting function for the infinitely long 



strip graph of type Gg, W{{Gs},q) defined in eq. (1.1), and to characterize its analytic properties in the 
complex q plane. We do not connect the ends of the strip in the longitudinal direction, i.e., the strip has free 
boundary conditions in both the transverse and longitudinal directions. As in Refs. pl|-p^, we let B denote 
the continuous locus of points in the q plane where W{{Gs},q) is nonanalytic. For the infinite-length limit 



of the strip graphs Gg defined in eq. (1.2) this locus of points B consists of arcs which do not separate the 
complex q plane into different regions. We calculate these for a number of different types of strips. 

In general, in the limit as the number n of vertices on an n-vertex graph of type G goes to infinity, the 
locus B forms as a merging together of zeros of the chromatic polynomial P{G, q) (called chromatic zeros of 
G). It is thus of interest to calculate chromatic zeros for strip graphs {Gs)m of finite length and to compare 
these with the asymptotic limit represented by the curves making up B. We perform this comparison for a 
number of types of strip graphs here. 

Our results are useful for several reasons. First, they are a set of exact calculations of chromatic polynomials 
for a large variety of strips, of both finite and infinite length. Second, the generating function introduced here 
enables us to calculate, in an algorithmic manner, the chromatic polynomials for very long finite strips. To 
appreciate how valuable such a method is, we recall the well-known fact |l^ that, in general, the calculation 
of the chromatic polynomial of an arbitrary graph is an NP-complete problem (NP = nonpolynomial) , i.e., 
roughly speaking, it requires a time which grows exponentially as a function of the number of vertices in 
the graph. Third, the generating function provides an elegant way to obtain directly the infinite-length 
limiting function Ty({Gs}, q) for a given strip graph of type Gg- These VF({G}, q) functions show a rich and 
fascinating analytic structure in the complex q plane. Fourth, as noted before, by studying this structure 
as a function of the width of the strips, one gains insight into how the limit of the infinite 2D lattice is 
approached (with free boundary conditions). 

We define 
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W[{G},q) = W^({G},g,)z,,„ = lim lim P(G,g)i/" (1.3) 

q — 'qs n — >oo 



and recall from Ref. pl[ | that at certain special points q^, the following limits do not commute: 



lim lim P(G',g)i/" 7^ lim lim P(G,(?)i/" (1.4) 

n — ^oo q — q — ^qs n — ^oo 



Our present definition (1.3) is the same as in Ref. and has the advantage of maintaining the analyticity 



of W^({G}, q) at the special points q^. 

In Ref. 1^ two of us had also discussed a second subtlety in the definition of M^({G}, q): for certain ranges 
of real q, P{G, q) can be negative, and, of course, when q is complex, so is P{G, q) in general. In these cases 
it may not be obvious, a priori, which of the n roots 

P(G,g)i/" = {|P(G,g)|i/"e2W"} , r = 0, 1, n - 1 (1.5) 

to choose in eq. ( |1 . l\ j . Consider the function VF({G}, q) defined via eq. ( pTl] ) starting with q on the positive 
real axis where P(G, q) > 0, and consider the maximal region in the complex q plane that can be reached 
by analytic continuation of this function. We denote this region as Pi. Clearly, the phase choice in ( |l.5| ) 
for g e Pi is that given by r = 0, namely P(G,(7)^/" = |P(G, (?)|^/". For the various exactly solved cases in 
Ref. concerned with families of graphs {G} for which the areas of analyticity of VF({G}, q) include other 



regions not analytically connected to Pi, there is not, in general, any canonical choice of phase in (1.5). 
However, for the strip graphs studied here, the curves B where W{{Gs},q) is nonanalytic do not separate 
the q plane into different regions, i.e., the entire q plane, except for the locus B, is in Pi, so this second 
point having to do with phase choice in (^^) is easily dealt with, as discussed below. 

This paper is organized as follows. In Section II we introduce the generating function. In Section III we 
calculate the generating functions for a number of interesting types of strip graphs Gs ■ Section IV is devoted 
to a study of the analytic properties of the resultant T/F({Gs}, q) functions. Some final remarks are given in 
the conclusions. Besides the mathematical references p"3|-p^ and the previous Refs. [|l]-|3 by two of us, 
some related work is in Refs. [|9[|6|j27|j2^||. We shaU report further resuhs in Ref. |9|. 



II. CALCULATION OF P{Gs,q) VIA A GENERATING FUNCTION 

A. Generalities 

It is possible to calculate the chromatic polynomials for a wide variety of strip graphs {Gs)m of the type 



defined in eq. (L2) by computing a generating function. This also enables us immediately to calculate the 
resultant asymptotic limiting function W{{Gs},q)- Before proceeding, we discuss some preliminaries. We 
denote the number of vertices in a graph G as n{G). Consider two successive repeating subgraph units H 
and denote their intersection as the graph 

Ln^Hr+inHr (2.1) 
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Let the number of vertices in the repeating subgraph unit H, the intersection graph Lh, and the initial 
graph I be n{H), n{LH), and n{I), respectively. Then the total number of vertices n in the strip graph 

{Gs)m is 

n{{Gs)m) = Cs,im + Csfi (2.2) 

where 

c.,0 = n{I) (2.3) 

and 

c«,i = n{H) - n{LH) (2.4) 
The values of these coefficients are listed in Tables I and II for the various strip graphs that we consider. 
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FIG. 1. Illustrations of strip graphs {Gi,(Ly))m of type s. Subscripts Ly and m refer, respectively, to the width and the 
number of repeated subgraph units H beyond the initial subgraph I (where I may be identical to H). For (a)-(g), m = 3, 
whereas m = 2 for (h). For ease of identification, the initial subgraph / is drawn with heavier lines. See text for further discussion 
of each type of strip graph, (a) (Gsq(3))3, an Lx x Ly = 5x3 strip of the square lattice; (b) (64^^(3) )3, an Lx X Ly = 5 X 3 
strip of the triangular lattice (represented as a square lattice with additional oblique bonds); (c) (G3q^(3))3, a strip of a square 
lattice along the diagonal direction, with the repeating subgraph H comprised of the upper, lower, and left members of a 3 X 3 
diagonally oriented box; (d) (G5q^(4))3, as in (c), but with a repeating unit H comprised of a diagonally oriented strip of 4 
squares; (e) (Ghc(3))3> ^ strip of the honeycomb lattice; (f) {Gkag{2))3y ^ strip of the kagome lattice; (g) (G4g8(2))3i ^ strip of 
the (4 ■ 8^) lattice; (h) (G4gg(3))2, a wider strip of the (4 ■ 8^) lattice along a direction rotated 45° relative to that in (g). 

For example, consider strips of the square, triangular, and honeycomb lattices, as shown in Figs. |l|(a,b,e). 
We represent the strip of the triangular lattice as a strip of squares with oblique bonds added, as indicated 
in Fig. |l|(b). For our present discussion we represent the honeycomb lattice of Fig. |^(e) as a brick lattice. 
Denote the length as vertices in the horizontal (longitudinal) direction, and the width as Ly vertices in 
the vertical (transverse) direction. This corresponds to a width of w — Ly — 1 layers of squares or hexagons 
and 2w triangles, for the respective lattices. The strips of the square and triangular lattices shown in Fig. 
|l|(a,b) have L^ = 5 and Ly = 3. These strips are thus L^ x Ly lattices with free boundary conditions in 
both X and y directions. The strip of the honeycomb lattice of Fig. ^e), when represented in brick form, 
has Lx — 9 for the top and bottom edges and L^ = 10 for the interior layer, with Ly = 3. Later we will 
study the limit in which L^ 00 with Ly fixed, i.e., strips of size 00 x Ly. The square strip has, as a 
repeating subgraph unit H , a vertical stack of w squares. Hence, starting at one end, which will be taken by 
convention to be on the right and to consist of an initial subgraph / which is identical to H, the elongation of 
the strip proceeds by the addition of successive vertical stacks on the growing (left) end. For the honeycomb 
and triangular lattices, respectively, H is an oblique stack of w hexagons and a vertical stack of 2w triangles, 
and the initial subgraph I — H . For the strips of the square lattice (along the row direction) and of the 
triangular lattice, the relation between L^ and m is thus L^ — m + 2, and 

n{H) = 2n{LH) = 2Ly for G.^^Ly) or Gt^L,) (2.5) 

while for the honeycomb lattice, 

n{H) - 2{2Ly - 1) , n{LH) - 2{Ly - 1) for Ghc(L,) (2.6) 

The number n(/) depends on one's convention for the initial subgraph; in these three types of strips, since we 
take I = H,we have n(/) — n{H) for each. The total number of vertices is ?^((Gs(i^))m) = L^Ly = {m+2)Ly 
for the square and triangular lattices, so 

Cs(L„),i = Ly , Cs(Ly),o = '^Ly for s = sq, t (2.7) 

while for the honeycomb (brick) lattice, n{{Ghc{Ly))m) — 2[(m + 2)Ly — 1], so 

Chc(L„),i = 2Ly , Chc(Ly)fi = 2{2Ly - 1) (2.8) 

It is straightforward to obtain the analogous formulas for other types of strip graphs G^; we list the results 
for Cs^i and Cs_o in Tables I and II. 
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B. Main Theorem and Proof 



The generating function for a strip of type Gs, denoted generically as T{Gs,q,x), yields the chromatic 
polynomials via the expansion 

oo 

r(G„g,a;) = ^P((G,)„,g)x" (2.9) 

ni—O 

where a; is a symbolic variable and P((Gs)m,(?) is the chromatic polynomial for the coloring, with q colors, 
of the strip {Gs)m with m repeated subgraph units H beyond the initial (m = 0) one, /. We find that the 
generating function T{Gs, q, x) is a rational function of the form 

T{Gs,q, x) = — — (2.10) 

V{Gs,q, x) 

with 

N{Gs.q,x)^Y.''-Mx' (2.11) 

and 

V{Gs.q.x) = l+Y,bs,k{q)x'' (2.12) 

k=l 

where the a^ i and hs^i are polynomials in q that depend on the type and width of the strip graph Gg- The 
polynomials as^i also depend on the subgraph forming the right-hand end of the strip (discussed further 
below), while the hs,i are independent of this subgraph. The respective degrees jmax and kmax of the 
numerator and denominator polynomials J\f{Gs, q, x) and ^{Gs, q, x), as functions of a;, depend on the type 
and width of the strip Gs- 

We can state these results formally: 



Theorem 1. 

Let a strip graph of type Gs be constructed by m successive additions of a repeating graphical subunit H 
to an initial subgraph / as in eq. ( |1.2| ). Then the chromatic polynomial for this strip graph can be expressed 
via the generating function (2.9) given by (^.10[)-(pTT2). 



Proof 

We use the deletion-contraction theorem, which states that the chromatic polynomial for a graph G 
containing two adjacent vertices v and v' connected by a bond (edge) is equal to the chromatic polynomial 
for a graph in which this bond is removed minus the chromatic polynomial for the graph in which these 
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vertices are identifiedQ Applying the deletion-contraction theorem to the initial subgraph /, we obtain a 
finite set of linear equations, with m-independent coefficients, for P((Gs)„i, q), P{{Gs)m-i, l), P{{Gs.e)m, l), 
and possibly P((Gs.e)m-i, <?), where Gs,e denotes a strip graph of type s but with a different end subgraph 
(e) at the starting point on the right. In general, the application of the deletion-contraction theorem yields 
several Gs,e's with different end subgraphs e. The set of equations is linear because the deletion-contraction 
theorem is a linear equation among chromatic polynomials. The set of linear equations is finite because 
the repeating subgraph unit H is a finite graph. Next, assuming the form ( ^.9|) for each type of strip and 
substituting into these equations, we obtain a resultant set of linear equations for the generating function of 
the original strip, r(Gs, q, x), and the set of generating functions for the strips with different ends produced 
by the action of the deletion-contraction theorem, T{Gs,e,q), which we solve. This completes the proof. □ 

We illustrate the method of proof by an explicit example. We first note that a graph in which each adjacent 
pair of repeating subgraph units H intersect each other in a complete graph^ Lh ~ Ki, has a very simple 
chromatic polynomial. Since the intersection of H with Lh is just Lh, the intersection theorem^ in graph 
theory implies that the chromatic polynomial P{H, q) factorizes according to 

P{H,q)^P{Ke,q)PriH,q) (2.13) 

where 

e-i 

P{Ke,q)^l[{q-t) (2.14) 

which defines the polynomial Pr{H,q). For a chain graph {Gs:Ki)m consisting of m repeated subgraphs H 
with an initial subgraph I = H, the chromatic polynomial is easily computed by the interative use of the 
intersection theorem: 

Pm,K,Uq) = ^p^^Tqr ^ (\i^^-^))Pr{H,qr^' (2.15) 
The generating function is 



The deletion-contraction theorem is an expression of the elementary fact the chromatic polynomial of G enumerates 
the coloring of the graph such that the colors of each pair of adjacent vertices are different. In particular, then, the 
number of colorings in which the colors on the pair v and v' are different is equal to the number of colorings in which 
the colors assigned to v and v' are allowed to vary freely (subject to the constraints from their other bonds) minus 
the number in which the colors assigned to v and v' are the same. Reading this equation in the opposite direction 
gives the addition-contraction theorem of graph theory, which we shall use later. See, e.g., Refs. |^|- ^] for further 
discussion of the mathematical background. 

^ The complete graph Kg is defined as the graph consisting of £ vertices, each of which is connected to all of the 
others with bonds. 

The intersection theorem states that if two graphs G and G' intersect in a complete graph Ki for some I, then 
the chromatic polynomial of the union of these graphs is given by P(G U G',q) — P(G,q)P{G' , q) / P{Ki, q). 
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For example, for a strip of m p-sided polygons (p-gons) with each adjacent pair of p-gons intersecting on a 
mutual edge (^^2), the chromatic polynomial is 

P((Gp_gon;Kjm, g) = qiq - l)Dp{qr (2.17) 

where 

and the chromatic polynomial for a circuit with p vertices is 

P{Cp,q) = iq-ir + i-inq-l) (2.19) 

SO that {q) ~ q — 2 and D4 (q) = q^ — 3q + 3. (In the following, for brevity, we shall generally write out the 
expression for D^{q) but not for Dk{q) for higher values of k.) Since the chromatic polynomial can easily be 
calculated for this type of strip without the use of our more powerful methods, and since, as will be discussed 
below, the points in the complex q plane where the resultant M^dGs}, q) function is nonanalytic are at most 
isolated branch points, so that B = $, i.e., the null set, we concentrate on more complicated types of strip 



graphs (e.g. of greater width) where the factorization in cq. (2.15) does not occur 
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FIG. 2. Graphs for square strips of width Lj, = 3 of type (a) Gsq,m=i (b) Gsq,e=st,m=i, (c) Gsq,e=tt,,n=l, (d) Gsq,e=t,m=i, 
and (e) Gsq,e=ttv.m—l- See text for discussion. 

Let us, then, illustrate our method of proof by considering a strip graph that consists of squares, Gsq- 
The method consists of using the deletion-contraction theorem on the end of the original strip graph, which 
generates also related strip graphs with different initial subgraphs; one then uses the deletion-contraction 
on these, thereby generating further types of strip graphs, and so forth until no further additional types of 
strip graphs are generated. (Of course, all of these additional types of strip graphs also have free boundary 
conditions.) In this manner, one obtains a set of linear equations for the generating functions of all of these 
types of strip graphs, which is then solved. A general strip of this type has size Lx x Ly vertices, i.e., 
(m + 1) X ui squares, where m — — 2 and w ^ Ly — 1, as given above. Both the strip with Ly = I 
(the line) and Ly — 2 have chromatic polynomials that factorize, so we consider the next simplest case 
of width Ly = 3 (see Fig. ||(a)). To simplify the notation as much as possible, we shall use the symbols 
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Ps{Ly,e),m = P[[Gs(Ly),{e))m,q) and sliall put Gsq{Ly=?,) = Gsq(3) = Gss- A first application of the deletion- 
contraction theorem on one of the vertical bonds (say, the upper one) at the right end of a length-m strip 
yields the equation 

Pss,.n = {q- l)Di{q)P,s,m-l - Pst,m (2.20) 

for TO > 1, where Pst.m refers to the chromatic polynomial for the graph Gsq^e=st,m with right-hand end 
consisting of a square and triangle (whence the subscript st) shown in Fig. |^(b). Given our convention of 
taking Gsq^m=o to be the graph with two squares on top of each other and Gsq,e=st,m=a to be a graph with 
a triangle adjoined to a square, it follows that 

Pss.o = P((Gs9(2))o, q) - qiq - 1)^4(9)' (2.21) 

and 

Pst,o = q{q-l){q-2)Diiq) (2.22) 



Equation (2.2C) thus relates the chromatic polynomial of the square strip of length to and width 2 to the 
chromatic polynomial of the same type of strip of length m — 1 and to the chromatic polynomial of this type of 
square strip with a different type of endgraph, as shown in Fig. ^(b). We next apply the deletion-contraction 
theorem to this latter type of strip, specifically to the remaining vertical end bond on the right in Fig. ^(b), 
to obtain the relation 

Pst^n. = {q- l){q - 2)Pss.m-l - Ptt,ra (2.23) 

for m > 1, and Pu.m is the chromatic polynomial for the square strip graph Gsq.e=u.m with a right-hand end 
consisting of two triangles, as shown in Fig. ^(c). Here we take Gsq,e=tt,m=o to be the two triangles adjoined 
along their horizontal edge, so that 

Ptt,o ^ q{q - l){q - 2f (2.24) 
Proceeding by applying the theorem to one of the hypotenuses of the triangles in Gsq,e=u,m, we derive 

Ptt,m — {q ^ '2)Pss,m-l — Pt,m (2.25) 

for TO > 1, where the graph Gsq,e=t.m=i is shown in Fig. |2| (d), so that the to = graph of this type is a 
single triangle, with Pt,Q = PiK^, q) = q{q — l){q — 2). The next deletion, on the round bond on the right 
in Fig. §(d), yields 

Pt,m — Pss,m—1 Pttv .m — 1 (2.26) 

for TO > 1, where the graph Gsq,e=uv,m=i is shown in Fig. ^ (e), and the to = graph of this type is taken 
to be the 5-vertex one at the right, with 
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Pttvs) = q{q - 1)!?' - + lOq - 7) 

Deleting one of the oblique bonds at the right end in Fig. ||(e) gives the relation 

Pttv.m — Di{q)Pss,m-l ^ {Q ^ '^)Pt,m 



(2.27) 



(2.28) 



for m > 1. Evidently, these linear relations have now closed, i.e., no new types of strip graphs are introduced 
beyond the ss, st,tt,t,ttv that we have. Next, we substitute the assumed form ( p.9[ ) into these equations. 
Note that we do not have to specify jmax or kmax] these are determined automatically by our solution of 
the equations. This substitution yields five linear equations for the five generating functions for the different 
types of strips: the original one, F^s, together with the other four that have been generated by the operations 
of the deletion-contraction theorem. This set of equations can be written as the linear transformation 



where the vectors and v„ are defined as 



Tst ^ 
V Tttv J 



(2.29) 



(2.30) 



/ Pss,0 + Pst,0 

Pst.O + Ptt,a 

Ptt.Q + Pt,0 
Pt.O 

V Pttvfi + {q- 2)Pt^o / 



and, with this ordering of the basis, the matrix T is given by 



T = 



/ l-{q-l)D4q)x 


1 











-{q-l){q-2)x 


1 


1 





-{q-2)x 





1 


1 





— X 








1 


X 


\ -D4iq)x 








(9-2) 


1/ 



(2.31) 



(2.32) 



The solution to eq. (2.29), Up — T^^Vg, yields the desired generating function Tgs, together with generating 
functions for the four other square strips of width 2 and different subgraphs on the right-hand ends. In 
particular, all of the five generating functions in have the same denominator V given as the determinant 



V = det(r) 

which is a general relation; in the present case, 

2 



(2.33) 



(2.34) 



fe=i 
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where 



6sg(3)a--(9-2)(<z'-3g + 5) (2.35) 

and 

&sg(3).2 = iq- l){q^ - + 13q - 11) (2.36) 

Note the general property that 

det(r(a; = 0)) = 1 (2.37) 

whence 

V{x = 0) 1 (2.38) 



which was aheady incorporated into our form (2.10)-(2.12) with &o = 1- For aU of the five strips, we get 



jmax — 2, i.e., the respective generating functions have the form 

where, as above, the subscript (e) denotes the initial (end) subgraph and is imphcit for the original strip, 
since its definition specifies the initial subgraph /. For brevity of notation, we shall sometimes also suppress 
the (Ly) subscript, understanding it to be part of the type s. For the original square strip itself, the terms 
in the numerator polynomial M are 



^sq 



(3)^0 = 9(9-1)^4(9)' (2.40) 



asq(3)A = -9(9 - l)'!?' - 69' + 139 - 11) (2.41) 

From the results above, it is straightforward to obtain the numerators for the other four generating 
functions of the strips with different subgraphs at the right. For example, the numerator for Tu = 
^{Gsq{Ly=3,e=u),q,x) has 

asq{Ly=3,e=u),o = Ptt,o = q{q - 1)(9 - 2)' (2.42) 

asq(Ly=3,e=ttU = -9(9 - 1)'(9 - 2){q - 3) (2.43) 

In Appendix 1 we present two alternate ways of calculating the generating function T(Gs, q, x) that are 
complementary to the method given above. 
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C. Some Basic Properties of Generating Functions 



In this subsection we prove some basic properties of generating functions using elementary properties of 



chromatic polynomials of arbitrary graphs. First, we observe that the zeroth order term in eq. (2.£) is just 
the chromatic polynomial of the initial subgraph /: 

r(G„ g, 0) = AfiGs^q, 0) = a,,o = q) (2.44) 

Next, we have 
Theorem 2 

T{Gs,q,x)^0 for q = 0, !,...£ -1 if G^^Ki (2.45) 

whence 

Af{Gs,q,x) = for g = 0, if Gs 2 Ke (2.46) 

asj{q) =0 for g = 0, 1, ...i - 1 V j if G, 2 Ke (2.47) 
where Ki denotes the complete graph with £ vertices. 



Proof 

One requires at least £ colors to color a complete graph Kg subject to the condition that no adjacent vertices 
have the same color (all vertices are adjacent on a complete graph), i.e., P{Ke, q) = for q = 0,1, ...£ — 1. 
It follows that if G is an arbitrary graph that contains at least one complete graph Kg, then P{G, q) = for 
q = 0, 1, ^ — 1. In particular, this applies for the special case G — Gg- 

Since Ki is just a vertex, and K2 is just a bond, and since all Gg graphs have bonds (indeed have no 
disjoint vertices), a general consequence of Theorem 2 is that 

T{Gs,q,x)^Q for g 0, 1 (2.48) 

whence 

M{Gs,q,x)^Q for q^Q,l (2.49) 



as,j{q) = , for g = 0, 1 V j 



(2.50) 



The case ^ = 3 in eqs. ( 2.45 )-( 2.47 ), where the subgraph is a triangle, K^, applies to our strip graphs of the 
triangular, kagome, and (3 • 12^) lattices. 
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Another general result follows from the fact that an elementary upper bound on a chromatic polynomial is 
P{G^ q) < q" (for integer q) since the right-hand side is the number of ways that one can color the n-vertex 
graph G without any constraint. Indeed, this is what motivates the use of the 1/n'th power in the definition 



of W{{G},q) in eq. (1.2). Since P{G,q) is a polynomial, this bound also applies for q generalized from 
integer to real and complex values. It follows that W{{G},q) is bounded above by W{{G},q) < q for real 
q, and, more generally, iVFdG}, g)! < |g| for complex q. It is thus convenient to define, as we have in our 



earlier work [^l|-g3|, the reduced function 

Wr{{G},q)^q-'W{{G},q) (2.51) 

which has a finite limit as \q\ oo. In the case of regular lattices {G} = A, in the context of large-q Taylor 
expansions of reduced W functions [p0|j2l|,|2^ it is also useful to define a related reduced function 



W{A,q)= ,,^^'^:?L (2.52) 

where C is the coordination number of the lattice (in graph theory terminology, the degree A of a vertex on 
the lattice) . The large-q series expansion in terms of the variable 

y=-l^ (2.53) 

is then 

oo 

W{A,y) = l + Y,Wny'' (2.54) 

n=l 

(This assumes that W{A,y) is analytic at y = l/q = 0; for a recent discussion of this, see Refs. [pT| , ^ .) 
We shall discuss these large-g series below in connection with our exact solutions for T/F({Gs}, q) for various 
strip graphs. 

For strip graphs {Gs,kp)m where Gg is fc-partite (fcp), i.e., can be uniquely colored with k colors such that 
no two adjacent vertices have the same color, the generating function satisfies 

TiGs,kp,q^k,x) = (2.55) 
1 — X 

reflecting the fact that 

PiiGs,kpU,k)^kl (2.56) 

III. RESULTS FOR EXPLICIT GENERATING FUNCTIONS 

In this section we give results of our calculation of explicit generating functions for various types of strips. 
Since some of these, in particular, those associated with heteropolygonal Archimedean lattices, involve 
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repeating subgraph units H more complicated than just single polygons, it is useful to give formulas for 
the number of vertices as a function of the length of the strip, as measured by the number, m, of repeating 
subunits added to the initial subgraph / forming the end of the strip. We do this in Tables I and II. 
These tables also contain information on the degrees of the polynomials a^.j and bs_k in the numerator and 
denominator of the generating functions for various types of strip graphs. As was evident in our proof 
of Theorem 1 and will be discussed further in section IV (see Corollaries 3.1 and 3.2), the choice of the 
initial subgraph / affects only M and not T) in the generating function and consequently does not affect the 
asymptotic limiting function W{{Gs\,q)- 



TABLE I. Properties of generating functions for strip graphs (Gs)m as defined in eq. (1.2) of width Ly, comprised of m 
repetitions of a subgraph H, with an initial subgraph /, as discussed in the text, and having jmax = 1 and k 

max — oee 

the text for the definition of the diagonal width d. The degrees of agi^];^^-^ ,j and b^i^i^^-^j^, as polynomials in g, are denoted 
deg{a^(^l^^'j j) and rfe(;(fe5(^^) j,), and the subscript s{Ly) is omitted to save space. Rs is the discriminant defined in eq. |4.14 



^b.p. is the number of branch point singularities in the degeneracy equation (4.1!:). Narcs and Ns. cares denote the number of 
arcs and self-conjugate arcs, respectively, in the continuous locus of points B where W{{Gs}, q) is nonanalytic. 



Gs 


Cs,l 


Cs,o = deg{ao) 


deg{a\) 


sq, Ly ^3 


3 


6 


7 


sqa, d = 3 


5 


4 


5 


tri, Ly = 3 


3 


6 


7 


(3-6-3-6),Lj, = 2 
(3-122), L^^2 


5 


6 


7 


10 


12 


14 


(4.8^), L, = 2 


8 


8 


10 



deg{bi) 


deg{b2) 


deg{Rs) 


Nt.p. 


arcs 


Ns.c. 


3 


4 


6 


6 


3 


1 


5 


6 


10 


10 


5 


1 


3 


4 


6 


6 


3 


1 


5 


6 


10 


8 


4 





10 


12 


20 


18 


9 


1 


8 


10 


16 


16 


8 






TABLE II. Properties of generating functions for strip graphs {Gs)m with jmax > 1, kmax > 2. Notation is the same as 
in Table 1. 



s jmax kmax Cs,i Cs,o = deg{ao) deg{ai) deg{a2) deg{as) deg{bi) deg{b2) degib^) deg{b4.) 



sq, Ly=A 


2 


3 


4 


8 


10 


12 




4 


6 


8 




he, Ly — 3 


2 


3 


6 


10 


12 


10 




6 


8 


6 




sqd, d = 4 


2 


3 


6 


8 


10 


12 




6 


8 


10 




tri, Ly = 4 


3 


4 


4 


8 


10 


12 


14 


4 


6 


8 


10 


{A-8^)Ly=3 


2 


3 


12 


16 


20 


18 




12 


16 


18 
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We next list the polynomials a^j and 6s,fc that we have calculated for several strip graphs Gs- Others are 
listed in Appendix 2. 



A. Strip of the Square Lattice 



For strips of the square lattice of width Ly, we take the repeating subgraph unit to be a vertical stack of 
w ^ Ly — 1 squares, and the initial subgraph / on the right to be identical to this repeating unit, I = H. 
In illustrating the proof of Theorem 1 , we have already discussed in detail our calculation of the generating 
function for the square strip of width Ly — 3 and have listed the resultant polynomials a^gf^^yQ, 05^(3) i, 
^s(j(3),ii E^^^d fosg(3),2 as cqs. (2.40), ( p. 41 ), ( p. 35 ), and (2.36). For the square strip of width Ly — 4, we find 



jmax — 2, kmax 



3, and, (with the notation G 



sq{Ly=4) 



G. 



<?(4); 



asg(4),0 = 9(9 - 1)^4(9)^ 

asgi4),i - -qiq - l)'(g' - 4? + 5)(2g4 _ I5q^ + 45q^ _ eSg + 40) 

a.,(4),2 = q{q - 1)'(9' - 16/ + 112g6 - 449g^ + 

llSOg* - 1829^-'^ + 1858g2 - lG84g + 279) 



sq(4),l 



7q-^ - 23q^ + Alq - 33 



5sg(4),2 



2g^ - 23q^ + 116g^ - 329q^ + 553q^ - 517q + 207 



&.q(4),3 = -q^ + 16g^ - 112g*^ + 449g^ - 1130g^ + 1829g^ - 1858g^ + 1084g - 279 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



B. Strip of the Triangular Lattice 

We consider a strip of the triangular lattice of width Ly . The H and / subgraphs were defined above in 
section IL The strip graph of this type with m repeated H subgraph units is denoted {Gt(^i^^-^)m- The case 
Ly — 2 yields a chromatic polynomial that factorizes, as in eq. (2.15). For Ly = 3 we find jmax = 1 and 
kmax = 2 with 



at(3),o = q{q - 1)(9 - 2)" 

at(3),i --9(9-l)'(9-2)'(9-3) 
&t(3)a = V + 792-189 + 17 
bti3).2 = iq-2f{q-3) 

The results for the next wider strip, with Ly — 4, are somewhat lengthy and are given in Appendix 2 



(3.7) 
(3.8) 
(3.9) 
(3.10) 
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C. strip of the Honeycomb Lattice 



The minimal-width strip for the honeycomb (also called hexagonal) lattice with a nonfactorizing chromatic 
polynomial is shown in Fig. |^(e) and has a width Ly = 3 (represented as a brick lattice). In contrast to 
the strips of the square and triangular lattice with the same Ly — 3 width, the honeycomb strip yields a 



generating function with numerator and denominator functions that are of higher degrees in x, viz., jmax = 2 
and krnax — 3. We find 

a/ic(3),o = q{q - '^)D&{qf (3.11) 

a/.c(3),i = -q{q - - ISg** + 76g^ - 262(7^ + 589g5 

-gOlg"^ + 943(7^ - GSlg^ + 266q - 47) (3.12) 

aM3),2 = 9(9 - 1)'(9 - 2)2 (3.13) 

&ftc(3),i = -q^ + - 285"* + 56g^ - 71g2 + 58g - 26 (3.14) 

Kc(3),2 = {q- ^f{q^ - lOg^ + 43g'* - lQ2q^ + 144^^ _ i2Qq + 49) (3.15) 

&M3).3--('Z-l)'('Z-2)' (3.16) 



D. Strip of the Kagome Lattice 

It is also of interest to calculate generating functions for the chromatic polynomials of strips of lattices 
which involve more than just one regular polygons. The uniform lattices that involve tilings of the plane by 
one or more regular polygons are called Archimedean (for a review, see Ref. [0); for some related studies, 
see Refs. lUIIQ) and, since each vertex is equivalent to every other, these lattices may be defined by the 



sequence of polygons that one traverses when making a circuit (say, in a counterclockwise direction) around 
any vertex, 



k = \{pT (3.17) 



where pi label the polygon and > 1 is the number of times that a given polygon occurs contiguously in 
the product. The total number of times that a polygon of type pi occurs is the sum of its contiguous and 
noncontiguous appearances. The Archimedean lattices include three homopolygonal (also called monohedral) 
cases: square, triangular, and honeycomb. The remaining Archimedean lattices are heteropolygonal, i.e., 
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are composed of more than one type of regular polygon. For the strips of heteropolygonal lattices, it is 
necessary to fix a convention to describe the width. We shall choose the convention of taking the width to 
be the number of layers of the polygon with the greatest number of sides, and, to retain a correspondence 
with the homopolygonal lattices, we shall label the polynomials in T(Gs, q, x) hy w + I (which is Ly for 
the homopolygonal strips). For our present study we consider first a strip of the kagome lattice, denoted 



(3 • 6 • 3 • 6) in the standard mathematical notation ( 3.17 ). We take the strip to be of width w — 1 and oriented 



as shown in Fig. |^(f). The initial subgraph / on the right is a hexagon, and the repeating subgraph unit H 
is the graph comprised of a hexagon and the two triangles adjacent to it. We find jmax — 1, k^ax = 2 and 

afeag(2),0 = lil - ^)D6{q) (3.18) 

afeag(2),i =-'7('Z-l)'(9-2) (3.19) 
bkagi2).i - -{q - 2)(9" - 6<7'^ + iV _ I6q + IQ) (3.20) 

bkagi2).2 = {q-l)Hq-^f (3.21) 

We have also calculated generating functions for two other heteropolygonal Archimedean lattices, viz., the 
(3 • 12^) and (4 • 8^) lattices. The results for these are given in Appendix 2. We next proceed to study the 
limit in which the length of a strip goes to infinity. 

IV. ANALYTIC STRUCTURE OF W{{Gs},q) 

To calculate W({Gs}, q) from the generating function T(Gs, q, x), we prove a basic theorem. Let us write 
the denominator function V of the generating function r(G's,(7,x) in factorized form, in terms of its roots 

'^s,k^ k — 1, kmaxi 

ViGs,q, a;) - 1 + £ hAq)^'' - 6.,fc„_ n (a; - ^sAl)) (4-1) 

k=l k=l 

or equivalently 

V{Gs,q,x) = fj (1 - A,,fe((z)x) (4.2) 
fe=i 

where 

A.,fe(<z) = (4.3) 
(and the roots satisfy the relation &s,fc„„^ (^ 1)'^""'^ Ilfc^r ^s.k — !•) 
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Theorem 3 



Consider a strip graph Gs as defined in eq. ( |l.2| ), constructed by m successive additions of a repeating 
graphical subunit H to an initial subgraph /. For a given value of q, denote the function Xs.kil) of maximal 
magnitude in eq. (4.2) as Xs,max{q)- Then the n ^ oo asymptotic limiting function W{{Gs},q) as defined 



in eq. (1.1) with eq. (1.3) is given by 



W{{Gs},q) = (A,,„„,(g))i/^- 



(4.4) 



where the definition of Cg^i was given in eq. ( p.2|) and a general formula for it was presented in eq. (2.4). The 
phase to pick in evaluating eq. ( |l.5|) and eq. ( |4.4[ ) is chosen so that for sufficiently large real q, W{{Gs}, q) 
is real. 



Proof 



First, for the case kmax ~ 1, where the chromatic polynomial for the strip factorizes according to eq. 



(2.15), one recognizes 



A,(g)--fe,,i(<z)-P,(i/,g) 



(4.5) 



where Pr{H,q) was given in eq. ( 2.13| ), and the result (4.4) then follows immediately from the elementary 
expansion 

•^(^-'^'").=AA(G.,<z,x)f:A.(,)'".™ 

m=0 



T{Gs,q,x) 



1 - Xs{q)x 



(4.6) 



together with eqs. ( p.S| ) and (|1T|), whence 

W({Gs},q) - hm (A,(g))"/" = (A,(g))i/=-^ 



(4.7) 



where we have used eq. ( |2.2[ ). Note that, for the generic case where Cs,i > 2, W{{Gs}, q) has branch point 
singularities at the zeros (if any) of Xs{q)- We shall next discuss the case k^ax = 2 and comment on the 
straightforward generalization to kmax > 3. For brevity of notation, we shall suppress the q-dependence of 
the Xg^kiq)- For k„iax — 2, eq. ( 2.10| ) becomes 



T{Gs,q,x) 



Af{Gs,q,x) 



^ Af{Gs,q,x) 

{I - Xs^ix){l ~~ Xs^2x) (As,l - As,2) 

MiGs,q,x] 



Xs,l 



1 - As, IS 1 - As, 2a; 



(As,i — As, 2) _ 



A 



m=0 



m+1 \ _m 



(4.8) 



(which is, of course, symmetric under As,i <-*■ As. 2-) Then, from the definition (1.1) with ( [l.3| ), taking the 



1/n'th power as n ^ 00, and again using eq. (p. 2D, we obtain the result ( [4.4| ). It is straightforward to 
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generalize this to the case k^ax > 3. In passing, we note that if some fixed subset N^q > 2 of the As,fc 



were identicaUy equal, this would not affect the result (4.4). (This is similar to the situation where some 
number of the eigenvalues of a transfer matrix for a spin model are identically equal ]3^ .) This is trivial if 
the identical As,fc's are not maximal; if they are, it is also clear from the fact that 

lim (7Ve9(A..™a.(g))™+M ^ (KmaAq))'/'"' , (4.9) 

independent of Neq. (This possibility is mentioned only for completeness; we have not encountered any strip 



graph Gs where any Xs^kil) are identically equal.) The choice of phase, r = in eq. (1.5) is determined by 
the fact that for sufficiently large real q, \s,max{q) > 0. 
This completes the proof. □ 

We next state some corollaries to this theorem. 
Corollary 3.1 



Given the strip graph Gs as defined in eq. (|l.2|), the asymptotic function W{{Gs},q) as defined in ( 1.1 



with (1.3) is independent of the numerator function A/'(Gs, q, x). It follows also, a fortiori, that the continuous 



locus of points B where W{{Gs},q) is nonanalytic, is independent of N'{Gs,q,x). 



Corollary 3.2 



Given the strip graph Gs as defined in eq. (1.2), the asymptotic function W{{Gs},q) as defined in (1.1) 
with (1.3) is independent of the initial subgraph / and depends only on the infinitely repeated subgraph unit 
H. Furthermore, a fortiori, the continuous locus of points B where VF({Gs}, q) is nonanalytic is independent 
of / and depends only on H. 

In particular, for this second corollary, recall that our proof of Theorem 1 showed that V was independent 
of the initial subgraphs of the strip graphs generated by the operation of the contraction-deletion theorem, 
as was noted before eq. ( 2.33| ). It is interesting to see Corollaries 3.1 and 3.2 in operation by studying the 



chromatic zeros of strip graphs of a given type Gs, varying the type of initial subgraph / and the number 
of repeated subgraph units m. For small m, we find, as expected, that the chromatic zeros clearly depend 
on what / one uses, but as m gets larger and larger, this dependence is progressively reduced, until it 
disappears altogether as m ^ cxd. Moreover, as we shall discuss in later work for strip graphs of the 
form ^(0^1 the asymptotic m ^ oo limiting functions W{{Gs}, q) do, in general, depend on both of 
the end graphs / and J. 

For a strip graph {Gs)m composed of m repeated subgraph units H adjoined to an initial subgraph I — H 
and having the property that the chromatic polynomial factorizes as in ( 2.15| ), the asymptotic limiting 
function is given by 
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(4.10) 



where we have used eq. (2.4) for Cs^i- For the proof of this result, we observe that from the formula (2.16), 
one can immediately identify {q) ~ Pr {H, q) , and eq. ( 4.10| ) then follows from Theorem 3 and its special 



case eq. (^) for k^ax = 1- For example, for strips of the square and triangular lattices of width Ly = 2, 



using the £ = 2 special case of eq. ( [1.10| ) , we obtain 

W^({G,,(2)},9) = (<z'-3g + 3)i/2 

and 

W{{Gt,2)},q) = q-2 



(4.11) 



(4.12) 



For the case k^ax ~ 2, the roots of ^{Gs, q, x) = as an equation in x are rs,fc, fc = 1,2, given by 
rs,± = i-bs,i ± V^)/(26s,2), so that the Xs,k, k = 1,2, are 



(4.13) 



where the discriminant Rs is 



R, = bl^- 46,,2 



(4.14) 



Corollary 3.3 



Given the strip graph Gg as defined in eq. (L2); the maximal Xs.maxiq) in eqs. (42) and (44) has the 
behavior, for large real q. 



Xs,max{q) ^ q""-^ as q 



(4.15) 



Proof 

A chromatic polynomial for an arbitrary n-vertex graph G is of degree n as a polynomial in q (with 



highest-degree term is g"). Hence, from the definition (1.1), it follows that W^({G}, g) ^ q as q ^ cx) for real 



q (more generally, |VF({G},g)| ^ \q\ for complex q). Combining this with Theorem 3, eq. (4.4) yields the 
result (|t|). 



We next prove a theorem which explicitly determines the continuous locus of points B where Vl^({Gs}, q) 
is nonanalytic: 

Theorem 4 
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Let a strip graph Gs be defined as in cq. ( |1.2| ). Then the continuous locus of points B where the asymptotic 



function W{{Gs},q) is nonanalytic is comprised of the points where two (nonzero) Xs.k terms in eq. (4.2) 
both become maximal and degenerate in magnitude, i.e., the solutions to the equation 

Xs,maxiq) = >^s,max'{q) (4.16) 

Proof 

This follows from the general analysis of W{{G},q), specifically part (d) of Theorem 1, given in Ref. 



Note that, from eq. (4.4), if Cs,i > 2, then W{{Gs},q) has isolated branch point singularities where 
s,max{q) = 0, but, in general, these do not lie on the continuous locus B. □. 

A simple corollary to this theorem is that if the strip graph Gs defined by eq. (|1.2|) has the property that 



its chromatic polynomial factorizes in the form of eq. (2.15), then M^dGs}, q), which we have calculated in 



eq. (4.1C), is analytic in the complex q plane, except for isolated branch points, so that B is the null set. 



0. This follows because in this case there is only one As, (recall eq. (4.6)) and hence no possibility for a 
degeneracy of maximal As./t's or resultant locus B. For example, in the case where if is a p-sided polygon 
and each successive p-gon intersects the next on a mutual edge so that Lh = K2, we have 

W({G„,;kJ, q) = (Dp(g))i/(p-2) (4.17) 



where Dp{q) was defined in eq. ( ^.18 ). The continuous locus B = % for these graphs. For p > 4, the function 



W {{Gp-gons:K2\ 1 q) bas isolated branch point singularities at the p ~ 2 zeros of the polynomial Dp(q), but 
these are not connected with any continuous locus of singularities B. Note that Dp{q) is a polynomial of 
degree p — 2; for odd p > 3, Dp{q) has a single real zero at q — 2 and {p — l)/2 conjugate pairs of complex 
zeros, while for even q, Dp{q) has only complex zeros. 

Before proceeding to discuss iy({Gs}, q) for specific strip graphs Gs, we recall that the elementary property 
that for any family of graphs {G}, the associated locus of points B is invariant under complex conjugation: 

B{q) = B{q*) (4.18) 

This follows from the fact that the coefficients of each power of 9 in a chromatic polynomial P(G, q) are real 
(indeed, integers) and the fact that B originates as a merging of zeros of P{G, q) as the number of vertices 

71 00 

V. W[{Gs},q) FOR SPECIFIC STRIP GRAPHS 

We next apply our general theorems to calculate the exact W{{Gs}:q) functions for specific infinitely 
long strip graphs. The continuous locus of points B where W^dGs},^) is nonanalytic is shown in each of 
the accompanying figures. As was discussed above, as the length of the strip graph goes to infinity, the 
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locus B forms by the merging together of the zeros of the chromatic polynomial [i.e., the chromatic zeros) 
of this strip graph. One gains insight into how this limit is approached by calculating chromatic zeros for 
reasonably long finite strips of each type. Accordingly, we show these in the figures. There are also isolated 
chromatic zeros that occur at q = = 0, 1, and, for graphs that contain one or more triangles, at Qs = 2. 
As discussed in Ref. |^ and the Introduction, in cases (i) where these occur via factors that scale like the 
lattice size, i.e., {q — gs)*", ^ > li they remain as zeros of the limiting function W{{G}, q); (ii) if, on the 
other hand, they occur as factors like (g — g^)" where u does not scale like n, then they do not appear in 



W{{G},q). (If one were to use the alternate definition based on the opposite order of limits in (1.4), then 
these zeros at q — qs would appear as discontinuities in W{{Gs}, q).) In case (ii) the different orders of limits 



in eq. (1.4) do not commute, whereas in case (i) they do. Depending on the type of strip graph, isolated real 
chromatic zeros may occur at certain other points also. For example, in Fig. |^(b) there is a chromatic zero 
near to, but not precisely at, g = 2 and another at g ~ 2.239, while in Fig. ^(b) there is a chromatic zero 
at g = (3 + V5)/2 — 2.6180..., etc. The simplest types of strip graphs to consider are those with chromatic 
polynomials that factorize in the manner of eq. ( 2.15 ), corresponding to kmax = 1; we have already discussed 
these above. 
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FIG. 3. Analytic structure of the function W({Gsg(^^) }, q) for Ly = (a) 3, (b) 4, where {Gsq(^^)} denotes the 
{Lx = oo) X Ly strip of the square lattice. The arcs are the continuous locus of points B where this function is nonana- 
lytic. For comparison, the zeros of the chromatic polynomial P{{Gsq^Ly))m,q) for (a) Ly = 3, m = 16 (hence n = 54 vertices) 
and (b) Ly = A, m = W (hence n = 72) are shown. In this and subsequent figures, the initial subgraph / is specified in the 
text. 



We proceed to consider strips graphs with denominator functions with kmax = 2, i.e., ^{x) that are 



quadratic in x (hsted in Table 1). For these, we have the exphcit results (4.13), (4.14). We first consider 
the strip graph of the square lattice of size x Lj^ = oo x 3. (As was noted at the beginning, we always 
use free boundary conditions in this paper.) From eqs. ( ^.35 ) and ( ^.36 ), we calculate the W function by 
starting on the real axis with q > 2 and analytically continuing from there. In this range, &s(;(3),i < smd so 
A. 



sg(3), 



Asg(3),-i where As,± was given in eq. (4.13). Hence, in this range 



W^({G.,(3)},g) = 2-1/3 



{q-2){q^-3q + 5) + 



{q^ -5q + 7){q^ - 5q^ + llq^ - 12q 



1/21 1/3 



(5.1) 



In Fig. |^(a) we show our calculation of the exact continuous locus of points B where W {{G sq(3)} , q) is 
nonanalytic. These consist of two complex-conjugate (c.c.) arcs together with a self-conjugate arc that 
crosses the real axis at g = 2, where &sg(3)^i vanishes. Since B does not enclose any regions, we can 



analytically continue the expression (5.1) throughout the q plane, exclusive of the points on the locus B 
itself. The two complex zeros of 6sg(3),i at g = 3/2 ± (i/2)VTT) = 1.5 ± 1.66z lie on the upper and lower arcs. 



The discriminant of the square root in eq. (5.1) has zeros a.t q = 5/2 ± {i/2)^/3 (for the quadratic factor) and 
q = 0.587 ± 1.14i and 1.91 ± l.lOi; (for the quartic factor). These zeros are all simple and hence are branch 
points of the square root. As is evident in Fig. ^(a), the first pair form the endpoints of the self-conjugate 
arc crossing the real axis, while the second two pairs form the endpoints of the two complex-conjugate arcs. 
In Fig. ||(a) we also show our calculations of the chromatic zeros for a long finite strip of the square lattice^. 
As expected, aside from the generally occurring isolated zeros at q — and 1, these zeros lie close to the 
arcs that comprise B in the infinite- length limit. 

Next, in Fig. ^(b) we present our exact calculation for B for a square strip of the next greater width, i.e., 
size X Ly — oo X 4:, together with chromatic zeros computed for a long finite strip. Our explicit analytic 
expression for W{{Gsq{4:)}, q), which involves the root of a cubic equation, is rather complicated and hence is 
not given here; it can be obtained from eqs. (3.4)-(^!^). The comparison of the results presented in Fig. ||(b) 
with those for the Ly = 3 strip in Fig. |^(a) shows how this locus B changes as the width of a strip of a given 



^ We recall the well-known fact that calculations of the zeros of chromatic polynomials are difficult because for a 
graph G with n vertices, the coefficients pi in the chromatic polynomial P{G,q) = X^"_]^Pn9" range from pri = 1 to 
much larger values for coefficients of terms with I ~ n/2, as is clear from the theorem that \pr\ > ("Zi)- 
this reason, it is necessary to use a large number of digits of accuracy in routines to solve for the roots. 
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type increases. We see that all three of the arcs elongate, and the single self-conjugate arc of the Ly — 3 strip 
develops a more complicated structure for the Ly ~ 4 strip. In particular, multiple points appear on this arc 
for Ly — 4. Here, we use the term "multiple point" in the technical sense of algebraic geometry, meaning a 



point on an algebraic curve where two or more branches of the curve meet 1 34 3^ . The index of the multiple 
point is the number of branches that meet; in Fig. ^(b) there is (i) a multiple point at g ~ 2.265 where a 
vertical component of B intersects a small line segment in B on the real axis, and (ii) two multiple points 
forming a complex-conjugate pair, at which three components of B meet. It is also interesting to compare 
the exact arcs that we have calculated for infinite strips, and the associated chromatic zeros for finite strips, 
with the chromatic zeros calculated in Ref. |^ for an La, x Lj, = 8 x 8 square lattice with the same (free) 
boundary conditions. Already for a width Ly — 4, the chromatic zeros for the strip resemble those of the 
larger patch of the 2D square lattice. There are branchings in the chromatic zeros for the 8x8 lattice at 
complex-conjugate points close to the locations where we find multiple points on B for the cx) x 4 strip. Our 
previous results on W{{G},q) in Ref. |2|] have shown that boundary conditions can strongly affect B and 
the associated chromatic zeros that merge to form this locus in the n — s- oo limit of an n-vertex graph. For 
example, B = $ for the infinite limit of a linear graph with free boundary conditions, whereas for the same 
graph with periodic boundary conditions, i.e., the infinite-n limit of the circuit graph, we found that B is 
the circle |g — 1| = 1, which separates the complex q plane into two distinct regions^. Nevertheless, it is 
intriguing to compare our results with the chromatic zeros computed for an 8 x 8 lattice with cylindrical 
boundary conditions, i.e., periodic in one direction and free in the other po| . We observe that there are 
cusp-like structures in the chromatic zeros on the 8x8 lattice in the complex-conjugate regions where there 
are gaps between the arcs on the oo x 4 strip. 

Using our results for W^({Gsg(i^)}, g), we can calculate the large-g series for the reduced fimctions 
^r{{Gsq(Ly)}^'l) defined for a general graph G in eq. ( 2.51 ). In order to compare with the large-q se- 
ries expansion for the infinite square lattice, we actually consider the function W{{Gsq(Ly)}i q) defined in 
eq. d^H). (For this purpose, we formally substitute C = 4, although not all of the vertices of the strip 
graphs have the same degree.) With y — l/{q — 1) as in eq. ( 2.53 ), we find W {{G ^^(^i)} , q) = 1 + y, 
Wi{Gsgi2)}, q) = il + y){l-y + y2)i/2^ whence 

W({G,,(2)}, ,) = 1 + iy - + + |/ - ^y - Hly^ - '^y^ + 0(y«) 

~ 1 + O.by - 0.125y2 + 0.5625y^ + 0.21092/'' 



® In contrast, in statistical mechanics, the existence of the thermodynamic limit entails independence of boundary 
conditions, except in the well-understood sense that they can determine which of degenerate long range orderings the 
system takes in a low-temperature, broken-symmetry phase. This is one of several differences between the behavior 
of W{A,q) and the somewhat analogous function for a spin model in statistical mechanics, exp{f{A,K)), wher e / i s 
the per site free energy on a lattice A, and K — (3J as in the text. Another is the noncommutativity of eq. (1.4), 
which has no analogue in statistical mechanics. See Ref. |2l| for further discussion. 
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-0.0352?/^ - 0.1143y^ - 0.06592?/^ + 0{y^) (5.2) 

and 

T77/r^ N . 1 1 2 59 3 44 4 32 , 613 6 5666 , ^, g, 

T^({Gs,(3)}, 9) = 1 + - + + ^y' - ^y' - + -^y' + 0{y') 

~ 1 + 0.333?/ - O.llly^ + 0.728y^ + O.lSly"* 

-0.044y5 - 0.093/ + 0.288y^ + 0{y^) (5.3) 

These may be compared with the first few terms of the series for the infinite lattice |]30t| 

W{sq,q) = l + y^ + y'^ + 0{y^) (5.4) 

We see that as the width of the strip increases, the coefficients in the small-y series expansion of the 
corresponding W function approach those of the infinite lattice, up to terms of order 0(jjP), where p ^ Ly. 
In particular, the coefficients of the terms y'^ , k — 1,2, 4, 5, 6 are smaller than those of the y^ and y'' terms. 
Evidently, as Ly oo, these smaller coefficients eventually decrease to zero, while the coefficients of the y^ 
and terms increase to their infinite-lattice values of unity. For progressively larger strips, similar behavior 
occurs for the coefficients of higher terms in the small-j/ series. In this way one can see from small-y series 
expansions of W functions on infinitely long strips how the infinite-lattice series arise. 

It is also of interest to see how the values of the W{{Gsq(Ly)}, o) functions for the strips of various widths 
Ly compare with those of the infinite square lattice, W {{G sq(^i^^^ao) \ ^ q) = ^{sQ, q) for various values of q. 
For example, for q = 3, we have W^({G,,(i)}, 3) = 2, W^({G,,(2)}, 3) = ^/3 = 1.73205..., iy({G,,(3)}, 3) = 
[{l/2){5 + VT7)]^/^ = 1.65846..., and I^({G,,(4)}, 3) = 1.624945... These decrease monotonically toward the 
infinite-lattice value ^ M^({G^,(oo)}, 3) = W{sq,3) = (4/3)^/2 ^ 1.5396O.... Defining the ratio 

W{{Gsqioo)},q) 

we have Rw{{Gsqii)},3) - 1.2990..., Rw{{Gsgi2)},5) - 1.12500..., Rw{{Gsqi3)},5) = 1.07720..., and 
Rw {{G sq(4)} , = 1.05543.... Thus the infinitely long strip graph of the square lattice with width Ly = 4 
already yields a W function which is, for q = 3, rather close to the value for the infinite square lattice. We 
have verified that this is also true for higher values of q, using our Monte Carlo measurements of W{sq,q) 
in (Table 1 of) Ref. We have also verified that strips graphs of other lattices show similar behav- 

ior. We have found that, for a fixed q > qd-^), Rw{{Gs(Ly)}, o) decreases toward 1 as the strip width 
increases. Here, following our discussion in Refs. plp^p^ , gc(A) denotes the maximal (finite) real value 
of q where W{A,q) is nonanalytic; in the above works, we obtained the values qdsq) — 3, qdt) = 4, and 
qc{hc) = (3 -I- \/5)/2 = 2.618.. for the square, triangular, and honeycomb lattices, respectively (see further 
discussion below on boundary conditions). 



Rw{{Gs,iLy)}. q) = \ (5-5) 
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FIG. 4. Analytic structure of the function W({G5q^(i^)}, g) for = (a) 3, (b) 4, where {Ggq^^^^)} denotes an infinitely 
long strip along a diagonal direction of the square lattice, for the case where the repeating subgraph H has diagonal width 
d = Ld- Zeros of the chromatic polynomial P{{Gi,q^(La))'''^' l) (^) -^d = 3, m = 9 (so n = 49) and (b) Ld = 4, m = 8 (so 
n = 56) are shown. 

In Figs. 1(a) and 1(b) we show our exact calculations of B for infinitely long strips of the square lattice along 
the diagonal direction. Finite strips of these two respective types were shown in Fig. |^(c,d). Comparing the 
two different-width strips, one sees that the two pairs of complex-conjugate arcs in the narrower strip merge 
together to form a single such pair on the wider strip, while the self-conjugate arc elongates. Furthermore, 
from a comparison of these two figures with the previous two, one sees several similarities but also observes 
that even for a given lattice, the details of B depend, in general, on the orientation of the strip with respect 
to the principal lattice directions. 
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FIG. 5. Analytic structure of the function W({Gt(^^)}, q) for Ly = (a) 3, (b) 4, where {G 
strip of the triangular lattice. For comparison, the zeros of the chromatic polynomials P{(Gf(^i^^-'i)m,<i) for (a) L 
(hence, n = 42) and (b) Ly = 4, m = 12 (hence n = 56) are shown. 



t(-Ly)} denotes the {Lx = oo X Ly 
3, m = 12 



In Figs. |(a) and |(b) we show similar results for strips of width Ly = 3 and Lj, = 4 of the triangular 
lattice. One sees a significant change in B as the width increases from 3 to 4; the self-conjugate arc crossing 
the real axis at g = 2.56984... for Ly — 3 disappears, while what were previously the right-hand ends of 
the complex-conjugate arcs become multiple points of index 3, each sprouting two additional arcs pointing 
in toward the real axis. The patterns of chromatic zeros for the long finite strips also change in the same 
manner, so that, as expected, for each width, these zeros lie close to the respective asymptotic loci B. Our 
explicit result for W for the Ly ~ 3 strip is 



- 7q^ + 18q- 17 
q^ ~ Uq^ + 81q^ - 250q^ - 



U2q^ - 436q + 193 



1/2- 



1/3 



(5.6) 



for real q to the right of the crossing at q ~ 2.56984 noted above, with appropriate analytic continuation to 
the rest of the q plane. 
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FIG. 6. Analytic structure of the function VF({Gftc(3) }i 9): where {Gfi(.{3)} the Lx x Ly = oo X 3 strip of the honeycomb 
(brick) lattice. Zeros of the chromatic polynomial P((Ghc(3))6i ?) {n- = 46) are shown. 

In Fig. I we show our calculations of S for a La; x oo x 3 strip of the honeycomb (brick) lattice. 

As was true with the square and triangular lattices, the arcs help one understand the pattern of chromatic 
zeros calculated for an 8 x 8 lattice with the same free boundary conditions in Rcf. the positions of 
these arcs are close to the positions of the chromatic zeros, both for long finite strips and the 8x8 patch. 
Furthermore, the cusp-like structures in the latter case occur approximately where (i) our exact arcs have 
multiple points and (ii) there are gaps between the arcs. 

Using our calculations of the generating functions, we have also determined the exact W functions and their 
associated continuous loci of nonanalyticities B for infinitely long strips of three heteropolygonal lattices. 



viz., kagome {— (3 • 6 • 3 • 6)), (3-12^), and (4 • 8'^), in the standard mathematical notation of eq. (B.17 ). As an 
example, we consider the kagome lattice. In Fig. |^ we show our exact determination of B for an infinite strip 
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of this type, a finite section of which was shown in Fig. [|(f). There is a particular feature of this calculation 
that is worth noting: for strip graphs which have T) functions that are quadratic in x, the As^fc((7)'s given 



by eq. ( 4.13 ) do not, in general, factorize. However, for the strip of the kagome lattice considered here, the 
As^fc((7)'s do factorize. This happens because fofeag(2),i a-nd the square root of the discriminant both have a 
common factor, (ci — 2). Related to this, not all of the zeros of the discriminant are branch point singularities 
of the square root of this discriminant; the zero at g = 2 is of even order (a double zero), and hence is 
not associated with a branch point of the square root of the discriminant (although it is a branch point 
singularity of the full W function; see below). We calculate 



Nl/5 



- 6q-^ + Uq^ - 16q + 10 



q^ - 12q^ + 6Aq^ - 200q^ + 404g^ - 548g^ + 500?^ - 292g + 92 



1/2- 



1/5 



(5.7) 



As before, one can obtain this for positive real q and analytically continue the expression to the full complex 
q plane, except for the locus B where W{{Gi;ag{2)}, q) is nonanalytic. Note that although the point q ^ 2 is 
a branch point singularity of W{{Gkag{2)}j 9)j it is an isolated singularity, and does not lie on the continuous 
locus of nonanalytic points B. 
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FIG. 7. Analytic structure of the function VF{{Gj.^(2) }i <?)> where {Gkag(2)} denotes an infinitely long strip of the kagome 
lattice, a finite-length section of which is shown in Fig. Mf). Zeros of the chromatic polynomial P{{Gkag(2))'.), <?) of ^ finite strip 
with m = 9 (hence n = 51) are shown. 



We present our results for the generating functions of the other heteropolygonal Archimedean lattices in 
Appendix 2 and the resultant determinations of B in the figures of Appendix 3. The chromatic zeros for 
long finite strips of each of these types are also shown in these figures. As with homopolygonal lattices, 
the arcs for the w = 1 {3 ■ 12^) and w ~ 2 {4 ■ 8^) infinite strips show close similarity in position not only 
with the chromatic zeros calculated for finite strips of the same type and width but also with our previous 
calculations of chromatic zeros for reasonably large patches of the respective lattices with Lx — Ly in Rcf. 
[Q. For the (3 • 12^) lattice, the large number (nine) of closely spaced arcs in Fig. [s] correlates nicely with 
the rather smooth and round distribution of chromatic zeros that we presented in Fig. 8 of Ref. For 
the (4 • 8^) lattice, the correspondence with the chromatic zeros of Fig. 9 in Rcf. is close for the w = 2 
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strip, shown in Fig. |9[ here, as before, one sees a faithful correspondence between cusp-hke structures in the 
chromatic zeros on the 2D path and the multiple points and gaps between arcs on B for the strip. 

There are a number of interesting features exhibited by our results. A general feature is that for all of 



these strip graphs of the form Gs ~ (YiiLiH)! as in eq. (1.2), the respective loci B form arcs which do 
not separate the complex q plane into different, disconnected regions. This finding may be compared with 
our exact results in Refs. |^,^. There we showed that the nature of the locus B depends strongly on the 
topology of a graph, (boundary conditions for a lattice graph) . In the simplest case of an n- vertex tree graph 
r„, in the limit n oo, VF({r}, q) is analytic in the entire finite q plane (i.e., B = the null set), whereas 
if one ties together the ends to make an n-vertex circuit graph C„, then in the limit n — )■ oo, VF({C}, q) is 
nonanalytic on the boundary B formed by the circle |<z — 1| = 1, which divides the q plane into two regions 
pl[ . Similarly, the n — > oo limit of an n- vertex strip graph (chain) composed of m p-sided polygons, denoted 
{Ch)p^n in Ref. (with n = {p — 2)m + 2) yields a W{{{Ch)p}, q) function with B ~ whereas if one 
ties together the ends to form a circular ladder graph, the resultant function W^({L}, q) is nonanalytic on a 
continuous locus B which forms the boundaries of four different regions in the q plane (as shown in Fig. 3 of 
Ref. [|2^). In all of these cases, B was cither (i) the null set 0, for the n oo limit of n- vertex graphs such as 
the complete graph Kn, tree graph T„, and strip (chain) graph composed of p-gons, (C/i)p^„, or (ii) a set of 
boundaries that separated the complex q plane into two or more different regions, within which W{{G},q) 
was analytic, such as for the circuit, C„, periodic and twisted ladder graphs, and others like the wheel graph 
and its generalization to the p-wheel graph pl| , p3t f[ Thus, the strip graphs studied here yield W{{Gs},q) 
functions with analytic structure intermediate between type (i) and (ii) found earlier: B is nonempty but 
does not divide the q plane into separate regions of analyticity. Moreover, different strips exhibit a wide 
variety in the number of arcs. As we shall discuss in Rcf. |2^], we have found that for strip graphs of the form 
"^(n"=i ^ contains arcs again, but these can actually enclose regions; this happens if the degeneracy 
equation ( [f.l6[ ) involves at least one Xs^k which is a polynomial in q. In contrast, for all of the strip graphs 
of the form (Jlfli ^^d having B ^ % considered in the present paper, the maximal As^^'s that occur in 



the degeneracy equation (4.16) are algebraic, but not simply polynomial, functions of g, i.e., they involve 
branch point singularities. 

The origin of this arc structure of B for the strip graphs studied here can be seen explicitly by examining 
the calculation of this locus in the simple case where k„iax — 2, i.e., where the denominator V^GsTqjx) 
of the generating function T(Gs, q, x) is quadratic as a polynomial in x. Here, aside from the trivial case 
where for some q, bs^i{q) — bs.2{q) = n the continuous locus of points B where W{{Gs},q) is nonanalytic 



The n-vertex wheel graph {Wh)„ is the graph formed by adjoining a point (i.e., connecting it with bonds) to all 
vertices of an (n— l)-vertex circuit graph C„-i. This may be denoted {Wh)„ = Ki x C„_i, where Kt is the complete 
graph on £ vertices. The generalized p- wheel graph is given by (Wh)^'' = Kp x Cn-p with {Wh)^'' = {Wh)n [ pl|j23| . 
® The restriction to nonvanishing bs^k is relevant for the T> functions for the kagome (i.e., (3 ■ 6 ■ 3 • 6)) and (3 • 12^) 
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is determined by the special case of eq. ( 4.16 ) with Xs.k{q), k = 1,2 given by eq. ( 4.13 ). The arcs comprising 
B extend between, and end on, zeros of the discriminant. It foUows that Narc — -^b.p./2, where Narc and 
A^fc.p. denote the number of arcs and branch points, respectively. In Table 1 we have listed these for the strip 
graphs with k^ax = 2 considered here. An analogous connection holds between the arc endpoints and branch 



points in the maximal- magnitude A^^fc's appearing in the degeneracy equation (4.16) for the case k^ax = 3, 
and one can see that it is a general feature of the degeneracy equations for As^^'s of maximal magnitude. 
Because B is invariant under complex conjugation, it follows that the arcs must occur as complex-conjugate 
pairs together with possible self-conjugate arcs. Clearly if Narc is odd, then there must exist at least one such 
self-conjugate arc. One can see by scanning Table 1 and the associated figures various examples illustrating 
both the presence and absence of self-conjugate arcs. 

The fact that B consists of arcs ending on branch points of the maximal-magnitude A^^fc in the complex 
q plane is reminiscent of the appearance of arcs (and line segments) of nonanalyticities of the free energy, 
ending on branch points of maximal-magnitude eigenvalues of the transfer matrix, in exact calculations of 



complex-temperature phase diagrams of ID systems (e.g., 35 3^]) and finite- width strips ||37|| . Our earlier 
work on complex-temperature phase diagrams |33|j3^ , |3^ helps to elucidate the conditions under which arcs 
or line segments do or do not enclose distinct regions in the respective complex plane (here, the q plane, there 
the complex plane in an appropriate Boltzmann weight variable). For example, in the ID nearest-neighbor 
Ising model with spin s = 1/2, the eigenvalues of the transfer matrix are given by \i/2,± — u~^^^(1±y^), and 
hence B consists of the semi-infinite line segment oo <u<0, where u = e'"^'^ with K = pj, f3 ^ l/iksT), 
and T and J are the temperature and spin-spin coupling. This locus of points ends on the branch point 
u — and does not separate the complex u plane into different regions. In contrast, for spin s = 1, the 
eigenvalues of the transfer matrix are Ai^i — u^^ ~ ui and 

Ai,± = (1/2) u^^ + l + ui± {u^^ - 2u{^ -I- 11 - 2mi -|- u\fl'^ (5.8) 

where u\ = e~^\ here the curves comprising B, which are the solutions to the degeneracy equation of 
maximal eigenvalues of the transfer matrix |Ai_i| — |Ai_+|, do separate the complex u\ plane into (four) 
different regions, as shown in Fig. 1(a) of Ref. [SHI. These curves end on the singular point at ui = and 



extend to complex infinity. Note that no arcs end on the finite branch points of the square root in eq. (5.S) 
because neither Ai_+ nor Ai,_ is a dominant eigenvalue at any of these points. Related to this, when one 
crosses one of the curves, one does not switch between two A's that are related only by changing the sign of 



strips of width Ly — 2; since both of these strips contain triangles, the chromatic polynomials vanish for all m > 1 
if g = 2. This is reflected in the property that for q = 2, the coefBcient functions as,i = = 6s, 2 = 0, and the 
respective generating functions become a;-independent. Hence, as is evident from Figs. Q and |^, the point g = 2 is 
not on the nonanalytic locus B. Related to this, as we have indicated in Table 1, in these two cases, although the 
respective degrees of Rs are 10 and 20, the number Nt.p. of branch points of \/R^ axe less by 2, viz., 8 and 18. This 
in contrast to the other strips listed in Table 1, where deg{Rs) = Nt.p.. 
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the square root, in contrast, e.g., to the situation described by eq. ( 4.16| ) for strip graphs considered here 



One also encounters arcs which close in such a manner as to separate the complex-temperature plane into 
different regions which are compact, rather than extending to infinity. This occurs, for example, in the case 
of the ID ^6 clock model ^] and the ID spin 1/2 Ising model with nearest- and next-nearest-neighbor 
spin-spin couplings Jnn and Jnnn for certain values of the ratio Jnnn / Jnn | |3^ . 

Another salient feature of our results is that although in some cases there is often a self-conjugate arc 
which crosses the real q axis farthest to the right, such an arc does not occur farthest to the left. This is also 
reflected in the reduced density (along arcs) of chromatic zeros near to q — (aside from the zero precisely 
at q — 0, which is always present for any chromatic polynomial). 

Our current results are relevant to the question of the conditions under which components of B are compact 
or, in contrast, extend infinitely far from the origin of the q plane, yielding a reduced function W{{G},q) 
which is nonanalytic at 1/q — Q. A family of graphs with noncompact B was studied in Refs. |p8| , pT| , and 
in Ref. infinitely many families of graphs were constructed with noncompact curves B, leading to the 
formulation of a general condition explaining when such noncompactness does or does not occur. Since for 
the Gs strip graphs considered here the branch points of the maximum-magnitude As,fc's, and hence the 
endpoints of the arcs, occur at finite values of q, these arcs never extend infinitely far from the origin q — 0. 
This is in accord with the criterion for compactness of B presented in our earlier general study ||2^ . 

A longstanding question in graph theory has concerned general features of chromatic zeros for a general 
graph G. Given the merging of a subset of these in the infinite-vertex limit to form the locus B of points 
where W{{G},q) is nonanalytic, one is also led to ask about general features of B for an arbitrary graph. 
We have discussed several of these above. Here we formulate a conjecture concerning the sign of the real 
part of chromatic zeros and thus of points on B. Before doing this, we recall that some time ago, it was 
conjectured [^8| that the chromatic zeros of an arbitrary graph have non-negative real parts. This conjecture 
was subsequently shown to be false, and now many examples of graphs are known with chromatic zeros 
having negative real parts. However, a striking feature of our exact calculations of chromatic zeros and 
nonanalytic loci B in Refs. [|l]j2^j2^ for many families of graphs, as well as the calculations of chromatic 
zeros for lattices with free, as opposed to cylindrical, boundary conditions, in Ref. [Q, is that they have 
Reiq) > 0. We are thus led to state the following conjectures. We first give a definition: a regular lattice 
graph has what we denote as a "global circuit" if and only if by travelling along one of the lattice vectors one 
can arrive back at one's starting point. The definition of lattice vector is clear for a cartesian d-dimensional 
lattice; for the triangular and honeycomb lattices, the lattice vectors can be taken to lie along 0°, 60°, and 
120°, etc. Thus, a patch of any of these lattices with free boundary conditions has no global circuits, while 
a patch with periodic boundary conditions in at least one direction has global circuits along this direction. 

Conjecture la 

If go is a chromatic zero of a regular lattice graph G with no global circuits, then i?e(qo) > 0. Furthermore, 
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the only chromatic zero with i?,e(go) = is go = itself. 
Conjecture lb 

Consider the n —f oo limit of a regular lattice graph G with no global circuits. The points q Cz B have the 
property that Re{q) > 0. Furthermore, the only point on B with Re{q) = is the point q = itself. 

Note that our exact calculations in Refs. |^,^ show that a graph may have global circuits and still 
have chromatic zeros and, in the infinite-n limit, a locus B, with the properties that (i) Re(q) > and (ii) 
Re{q) — Q — 0] this was true, in particular, of the circuit graph, its generalization [ p3| to Kp x C„_p. 
This is also true of the ladder graph (strip of the square lattice of width Ly = 2 with ends tied together (in 
an untwisted or twisted manner), i.e., with periodic boundary conditions in the longitudinal direction 
both of which have global circuits along the longitudinal direction. This shows that the condition of no global 
circuits is a sufficient, but not necessary, condition for properties (i) and (ii) above. Although our previous 
(exact) calculations of chromatic polynomials and W functions in Refs. showed that these are affected 

by boundary conditions, it also showed that, in all the cases considered, when one compared a given type 
of graph with different topology (i.e., existence or nonexistence of global circuits) due to different boundary 
conditions, there did exist a unique value qc such that for real q > qc, the W functions calculated with the 
different boundary conditions coincided. For example, for the line graph with free boundary conditions (and 
more generally the tree graph), in the infinite-vertex limit, iy({r}, q) = q—l, while for the circuit graph C„ 
{i.e. the line graph with periodic boundary conditions), in the same limit, iy({C}, q) = q — 1 for I? — 1| > 1 
and \W{{C},q)\ = 1 for |q - 1| < 1 (see section IV of Ref. |l|). Thus, W{{T},q) = W{{C},q) for real 
9 > 9c, where = 2 in this case (more generally, they also coincide throughout the region |g — 1| > 1 in 
the complex q plane). Similarly, for the n- vertex ladder graphs L„, in the n — > cx) limit with free boundary 
conditions, W{{L} fbc q) = {q^ — 3(7 + 3)^^^, while for (twisted or untwisted) periodic boundary conditions, 
this expression again holds in the region denoted Ri in Ref. | |2l[ |, defined in general as the maximal region 
which is analytically connected with the real interval q > qc for a given lattice^. For the periodic ladder 
graphs, qc is again equal to 2 (i?i also includes the interval g < on the real axis), and there are three 
other regions in the complex q plane where W takes on other analytic forms. These exact results suggest 
the following generalization: for a given type of strip graph Gg of a regular lattice, for real q > q^ where qc 
is the maximal (finite) real value at which W{{Gs},q) is nonanalytic (equivalently, B crosses the real axis) 
for some set of boundary conditions, the functions Vt^({Gs}, q) calculated with the different sets of boundary 
conditions coincide. A further generalization is that these W functions can be stated as follows. Consider 
the set of (in general different) loci B calculated for the infinite- vertex limit of a regular lattice with different 



® In eqs. (4.22)-(4.24) of Ref. the W function was defined per rung rather than per site; for the definition per 
site, as used here, the right-hand sides of these equations are raised to the power 1/2. 
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boundary conditions. Consider the maximal envelope of B for this set. Then the W functions agree in 
the region outside this maximal envelope, which region includes the area extending infinitely far away from 
the origin of the q plane, i.e., the inverse image of the origin of the 1/q plane. Note that this property is 
necessary for the existence of large-g Taylor series expansions of W for regular lattices. As discussed in Refs. 
1^,^, the existence of these large-g series expansions for W also depends on the property that in these 
cases B does not extend infinitely far away from the origin in the q plane. 

VI. CONCLUSIONS 

In conclusion, in this paper we have calculated the chromatic polynomials P{Gs, q) for strip graphs com- 
prised of repeated subgraphs H adjoined to an initial graph / and, from these, in the limit of infinitely long 
strips, the resultant W{{Gs}, q) functions. To do this, we have developed a powerful method in which we 
compute a generating function r(G's, q, x) which yields these chromatic polynomials as coefficients in a series 
expansion in x. We have studied the analytic structure of W{{Gs}, q) for a number of different types and 
widths of strip graphs. The locus of points where W{{Gs},q) is nonanalytic consists of arcs which do not 
separate the complex q plane into different regions. From comparisons of arcs for different-width strips of 
given types, we have shown how these arcs elongate and tend to move toward each other as the width of 
the infinitely long strips is increased. We have calculated chromatic zeros for long finite-length strips and 
have shown that, aside from the generally occuring real zeros at g = 0, 1 and, for strips containing triangles, 
(7 = 2, these lie close to the positions of the arcs. As we have discussed, these calculations of the ground 
state degeneracy W^(A, q) of the g-state Potts antiferromagnet on these infinitely long strips of various widths 
helps one to understand the behavior of this function on infinite 2D lattices A. In turn, this adds to our 
understanding of ground state degeneracy and entropy, a phenomenon which is exhibited in the physical 
world by such compounds as (II2O) ice. 

Acknowledgments 

This research was supported in part by the NSF grant PHY-93-09888. 

VII. APPENDIX 1 

In this Appendix we present two alternate ways to calculate the generating function T{Gs,q,x) that are 
complementary to the method given in section IIB. 
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A. Calculation of the Generating Function by Matching of Terms 



Here, as before, we solve for T{Gs, q, x) simply by matching the Taylor series expansion of T{Gs, q, x) 
to explicit calculations of the chromatic polynomials for the first few values of m. Since there are Np = 
jmax + kmax + 1 polynomials in (2.1C)-( 2.12| ), we can determine these by (i) calculating the chromatic 
polynomials P((Gs)m, q) for the various strips of type Gs with m ranging from to Np, (ii) performing the 
Taylor series expansion of r({Gs}, q, x) about a; = to order Np — 1, (iii) equating the successive terms in 



the Taylor series expansion to the respective chromatic polynomials, using (2.£), and (iv) solving for the aj, 
j = 0..,jmax and k = l..,kmax- For example, for many strip graphs of width w = Lj, — 1 = 2, we find 
that jmax — 1, k„iax = 2 0, SO that r{Gs,q,x) has the form ( 2.39 ) given above. Evaluating eq. (2J), one 
has 



Os.o = P{{Gs)o,q) 
(where P((G,)o,g) = P(/,g) from eq. ^A^) 

0-3,1 - asfibs,! = P{{Gs)i,q) 



(7.1) 



(7.2) 



and 



2as,obssbs,2 - assbs,2 + a^.i^s i - as.ob^ i = P{{Gs)3, q) 



which has the solution given by (7.1) together with 

2P((G,)o, g)P((G,)i, (7)F((G,)2, g) - P((G,)o, q?P{{Gsh, q) - P((G,)i, g)^ 



as,i 



bs,i 



Q 



P((G,)i, g)P((G,)2, g) - P((G,)o, g)F((G,)3, g) 
Q 



bs,2 = 



PiiGs)i,q)PiiGs)3,q) - PiiGs)2,q)' 



where 

Q = P((G,)o, g)P((G,)2, g) - P((G,)i, g)^ 
The resulting ai, bi, and &2 are polynomials, although this is not manifest in eqs. (^ 



(7.7) 



(7.3) 
(7.4) 

(7.5) 
(7.6) 
(7.7) 

(7.8) 



The honeycomb strip graph is an exception; for w = 2, the resultant generating function has jmax = 2, k max — O. 
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B. Calculation of the Generating Function by Addition-Contraction Theorem 



Another alternate method of calculation is based on the addition-contraction theorem from graph theory. 
This method is quite powerful and allows us also to calculate chromatic polynomials and asymptotic limiting 
functions W for strip graphs with specific subgraphs at both ends, not just one. We first recall the statement 
of the addition-contraction theorem: let G be a graph, and let v and v' be two non-adjacent vertices in G. 
Form (i) the graph Gadd. in which one adds a bond connecting v and v' , and (ii) the graph Gcontr. hi 
which one identifies v and v' . Then the chromatic polynomial for G is equal to the sum of the chromatic 
polynomial for the graphs Gadd. and Gcontr. ■ For our proof, with no loss of generality, it is convenient to take 



the initial subgraph / in eq. (1.2) to be identical to the repeating subgraph unit H; once having calculated 
the generating function for this case, it is straightforward to obtain the generating function for the strip 
with an initial subgraph / attached to the right end. By applying the addition-contraction theorem to the 
right-hand side of the initial H subgraph in (G(,)m, we obtain a set of strip graphs {G'g j)m with complete 
subgraphs labeled by j on the right-hand end. For example, let us consider a strip {G'g)m — TYiLi H with 
Ly = 4 vertices in the transverse direction, and label the vertices on the right-hand end, in sequence, as 
vitV2,v^,V4^. Now apply the addition-contraction theorem to the pair ui,W4; this yields the equation 

P((G'J„, q) = P((G;^,(i^4))„, q) + P((G;^,^^„J™, q) (7.9) 

where (G^ jjj--^ 4^)™ denotes the strip graph with the right-hand end modified by the addition of the bond 
connecting vertices vi and W4, and (G(,^^^„^)m denotes the strip graph with the right-end modified by 
identifying vertices vi and W4. Applying the theorem again to each of these two strip graphs, one obtains 
the equation 

P((G'J„, q) = P{{G',^Kjm,q) + P{{G',^K,iv^=v,))m,q) + P{{G',,K,iv,=.,))ra, q) 

+P((G ),n,q)+P{{G )rn,q) (7.10) 

Let us label the five resultant complete graphs on the right-hand end via the label i, with chromatic poly- 
nomials denoted by P{{G'g ^)o, q). Next, working one's way leftward from the right-hand end, apply the 
addition-contraction theorem in a similar manner to the next set of transverse vertices, and label the five 
resultant complete graphs with the label j. Define a (square) matrix M (in this case 5x5) with elements 
consisting of the chromatic polynomials of this set of graphs with complete subgraphs i on the right and j 
on the left. This procedure transforms the initial strip into a sum of factorized strips with parts that overlap 
in complete graphs. Therefore the intersection theorem yields, for m > 1, 

Pi{G'^^,U,q) = E ^'^''p((^^')tg)'''^^ = ^(Afi?),,P((G;j„_i,g), (7.11) 



where Z? is a (5 x 5) diagonal matrix with elements Di^i = 1/P((G(. i)a,q)- Because of ( 7.1C ), the generating 
function of Gl can be written as 
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r(G;,<7,a;) = ^r(G;^-,g,a:) 
j 

Writing the generating function of G'^ j in the more convenient form 

oo 
m— 1 

and using JtHI ) and d?!^ ), we obtain 

oo 

r{G'^,q,x) = Va(l - xMD)-^Mvb = va ^ a;"(AfL')"Afvb, 

m— 



(7.12) 



(7.13) 



(7.14) 



where in this case, for the strip G'^, Va = (1,1,1,1,1) and Vb — Va* denotes the transpose of Va. Note 
that in general Va and Vb can be chosen independently, and different choices of these vectors correspond 
to different boundary conditions on the two ends of the strip. This approach to the theorem allows one to 
write (|2.9|) in the equivalent form 



F((G'J,„,(7) =Va(Afi?)'"Mvb 



(7.15) 



The generating function (7.14) is a rational function of the form (2.1C) with the denominator given by 

ViG'„q,x)^l[il~Xr{q)x), (7.16) 

r 

where \r{q)'s are eigenvalues of the product of matrices MD. The boundary conditions on the two ends of 



the strip (defined by Va and Vb) determine which eigenvalues enter in the product in equation (7.16). Using 
these methods, we have also calculated chromatic polynomials and their asymptotic limits for various strip 
graphs of the form (Gs)m — J(YVeLi H)I (with (Gs)o = JI) where the strip has two, in general different, 
subgraphs on its ends. We shall present these calculations and related discussions of boundary conditions 
other than free ones in a subsequent paper p9t]. 



VIII. APPENDIX 2 



In this Appendix we list the a^j and bs^k polynomials making up the generating functions r(Gs,(7, x) for 
various strip graphs Gg. 



A. Strip of the Square Lattice in the Diagonal Direction 



Here we consider a strip of the square lattice such that the longitudinal axis of the strip is the diagonal 
direction of the lattice. In the simplest case, where adjacent squares touch each other only at a common 
vertex, the chromatic polynomial factorizes, as q{q— l)D4{q)"^^^ for m squares adjoined to an initial one, /. 
The first case with a nonfactorizing chromatic polynomial is for the strip shown in Fig. |l|(c) , where the initial 
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subgraph I at the right-hand end is a single diagonally oriented square, and the repeating subgraph unit H 
is a "wedge" consisting of two squares, one above the other, with a third to the immediate left. We denote 
the strip graph with m of these repeating units (G'sg^(3))„i. We find a generating function with jmax = 1, 



s<2d(3) 



aa9i(3),l =-g(g-l)''^('?-2) 

^ = -q5 _,_ _ 28g3 + 53^2 _ ^ 27 

^>s9,(3),2 = (g-l)(9-2)i?4(gf 



We have also calculated the generating function for a similar lattice of greater width, shown in Fig. |i|(d). 
For this case we find jmax — 2, k„iax — 3, and 



as9d(4),o = 9(9- 1)^4(9) 
asg,(4),l = -q{q - lfDi{q){3q'' - 25q^ + 86q^ - 152q^ + 136q - 47) 

a,g,i4),2^qiq-if{q-2)^D^{qf 

bsq^{i),i = -q^ + lOg^ - 459" + 1169^ - 185g^ + 178g - 82 
bsq,(i)a = D4{qf{3q^ - 22q^ + 6Aq^ - 88q + 49) 



8.1) 
8.2) 
8.3) 
8.4) 



8.5) 
8.6) 
8.7) 



8.9) 



fos9,(4).3--(9-2)'i?4(g)' (8.10) 

B. strip of the Triangular Lattice of Width Ly — 4 

In the main text, we have given our results for the generating function for a strip of the triangular lattice 
with width Ly = 3; here we list the results for the next wider strip of this type, having a width Ly — 4. We 
calculate that jmax = 3, kmax = 4, and 

at(4),o =«('Z-l)(9-2)6 (8.11) 
at(4),i = -q{q - l){q - 2f{q ~ 3)(3g'^ - 25q^ + 76g2 _ 98g + 43) (8.12) 
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at(4),2 = q{q - l)iq - 2)^(g - 3 f{3q^ - 12q + 8) (8.13) 

at(4),3 = -q{q - ifiq - ^fiq - 3)^ (8.14) 

&t(4)a = -q^ + Wq^ - 42g2 + 88q - 76 (8.15) 

&t(4).2 = iq- 2)(g - 3)2(3g^ - 22g2 + 60g - 60) (8.16) 

^t(4),3 = -{q - 2)'(g - 3)3(39^ - 21(?2 + 5l9 - 43) (8.17) 

&t(4).4 = (g-2)6(g_3)4 (8.18) 



C. strip of the (3- 12^) Lattice 

We consider a strip oriented such that the upper and lower sides of the 12-gons are parallel to the direction 
of the strip. The width of the strip is given by the diameter of one 12-gon and the initial subgraph / on the 
right is taken to be a 12-gon and the repeating subgraph unit H is a 12-gon with its two adjacent triangles 



(which do not touch each other). We calculate jmax — 1, k^ax = 2 and 

a3i2(2),o = 9(9 - 1)-Di2(g) (8.19) 

a3i2(2),i = -q{q - 1)'°(9 - 2)iq^ - 5g + 5) (8.20) 

^312(2), 1 = -{q - 2)(g^ - llg^ + 54g^ - 156g*^ + 294g^ - 378g^ + 336g^ - 204^^ + 82q - 22) (8.21) 

^3i2(2),2 = (9 - lYiq - 2)Hq^ - 5<7 + 5) (8.22) 



D. Strip of the (4 ■ 8^) Lattice 

The minimal width strip of the (4-8^) lattice is shown in Fig. |^(g). The initial subgraph / is an octagon 
and the repeating subgraph unit H is an octagon with two squares at the same x coordinate (which do not 
touch each other). For this strip we find jmax = 1 and k^ax — 2 and calculate 

a488(2),o = q{q - i)^8(g) (8.23) 

"488(2),! = -qiq - ifiq - 2)(3<72 -9q + 7) (8.24) 
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&488(2),i = -q^ + llg^ - 55q^ + 163q^ - 31 V + 406?^ - 354g2 + 203q - 63 (8.25) 



5488(2), 2 



= {q- l)^{q - 2)(3g2 - 9q + r)Di{qf (8.26) 



We have also calculated the generating function for the next wider strip of the (4-8^) lattice, shown in 
Fig. 0(h). The longitudinal direction of this strip is rotated 45° relative to that of the one above. We find 
jmax — 2, k„iax — 3, and calculate 

a488(3),o = qiq - ^)iq^^ - 18<7^^ + I53q^'^ - 814(7" + 3030?^° - 8358g^ + 17655g*^ - 29104/ 

+377965*^ - 38720g^ + 31057g^ - 19118?^ + 86669^ - 2649g + 427) (8.27) 

0488(3), 1 = -qiq - ^)*{q^^ - 23g" + 251g^^ - 1719g^^ + 8241^" - 29235g^° 
+79172g9 - 166558(7^ + 274466g^ - 354610(7^ + 356940g^ 

-275630g^ + 158575g^ - 64360^2 + 16492^ - 1999) (8.28) 

a488(3),2 = q{q - 1)"(<Z - 2f{q - ifDiiqf (8.29) 

&488(3),i = -9^^ + 16g" - 120gi" + 558g^ - 1794gS + 4212^^ - 7437g^ 

+ 10018g^ - 103249"* + 8064^^ - 4648^^ + 1854g - 414 (8.30) 

^'488(3),2 = {q- l)''(g^^ - 20g" + 188gi° - 1094g9 + 4375g^ - 12640^^ + 27033?*^ 

-43164g^ + 51235g* - 44380g^ + 26931?^ - 10462g + 2017) (8.31) 

&488(3),3 = -(9 - 1)'('7 - 2)2(r7 - ifDiiqf (8.32) 
IX. APPENDIX 3 

In this Appendix we include a number of additional figures depicting exact results on the analytic struc- 
ture for Ty({Gs},g) functions for heteropolygonal Archimedean lattices Gs, together with calculations of 
chromatic zeros for long finite strips of each type. 
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FIG. 8. Analytic structure of the function W({G(3.i22) }i 9)1 where {03.^22)} denotes an infinitely long strip of the (3-12^) 
lattice oriented as discussed in the text, and of width one layer of 12-gons. For comparison, the zeros of the chromatic polynomial 
-P('-'(3 i22))3i9) of a finite strip with m = 3 (hence n = 42 vertices), are shown. 
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FIG. 9. Analytic structure of the function W{{Gf^^,g2^f^2)}'l)±^'^^^^^ {^(4'8'^){w)} denotes an infinitely long strip of the 
(4 • S'^) lattice of width u>; a finite-length section is shown in Fig. bj(h). For comparison, the zeros of the chromatic polynomial 
P((G(-4.g2-)(3-))2i q) of a finite strip (with n = 40) is shown. 



46 



[1] W. F. Giauque and J. W. Stout, J. Am. Chem. Soc. 58 (1936) 1144). 

[2] L. Pauling, The Nature of the Chemical Bond (Cornell Univ. Press, Ithaca, 1960), p. 466. 

[3] E. H. Lieb and F. Y. Wu, in C. Domb and M. S. Green, eds.. Phase Transitions and Critical Phenomena 

(Academic Press, New York, 1972) v. 1, p. 331. 
[4] N. G. Parsonage and L. A. K. Staveley, Disorder in Crystals (Oxford, 1978). 
[5] P. W. Atkins, Physical Chemistry (Freeman, San Francisco, 1994), p. 710. 
[6] G. H. Wannier, Phys. Rev. 79 (1950) 357. 
[7] K. Kano and S. Naya, Prog. Theor. Phys. 10 (1953) 158. 

[8] A. P. Ramirez, G. P. Espinosa, and A. S. Cooper, Phys. Rev. Lett. 64 (1990) 2070; C. Broholm, G. Aepph, G. 

P. Espinosa, and A. S. Cooper, ibid. 65 (1990) 3173. 
[9] A. S. Wills, A. Harrison, S. A. M. Mentink, T. E. Ma son, and Z. Tun, "M agnetic Correlations in Deuteronium 

Jarosite, a Model S = 5/2 Kagome Antiferromagnet" , :ond-mat/960710C . 
[10] C. Zeng and V. Elser, Phys. Rev. B42 (1990) 8436; R. R. P. Singh and D. Huse, Phys. Rev. Lett. 68 (1992) 

1766; N. Elstner, R. R. P. Singh, and A. P. Young, ibid 71 (1993) 1629; N. Elstner and A. P. Young, Phys. Rev. 

B50 (1994) 6871. 
[11] R. B. Potts, Proc. Camb. Phil. Soc. 48 (1952) 106. 
[12] F. Y. Wu, Rev. Mod. Phys. 54 (1982) 235. 
[13] G. D. BirkhofT, Ann. of Math. 14 (1912) 42. 

[14] H. Whitney, Ann. of Math. 33 (1932) 688; Bull. Am. Math. Soc. 38 (1932) 572. 
[15] G. D. BirkhofT and D. C. Lewis, Trans. Am. Math. Soc. 60 (1946) 355. 

[16] R. C. Read, J. Combin. Theory 4 (1968) 52; R. C. Read and W. T. Tutte, "Chromatic Polynomials", in Selected 
Topics in Graph Theory, 3, eds. L. W. Beineke and R. J. Wilson (Academic Press, New York, 1988; R. Nelson 
and R. J. Wilson, eds.. Graph Colourings (Longman, Essex, 1990). 

[17] W. T. Tutte Graph Theory, vol. 21 of Encyclopedia of Mathematics and its Applications, ed. Rota, G. C. (Addison- 
Wesley, New York, 1984). 

[18] F. Harary, Graph Theory (Addison- Wesley, Reading, 1969); N. L. Biggs, Algebraic Graph Theory, (Cambridge, 

U.K., Cambridge Univ. Press, 1st ed. 1974, 2nd ed. 1993). 
[19] E. H. Lieb, Phys. Rev. 162 (1967) 162. 

[20] R. J. Baxter, J. Phys. A 20 (1987) 5241. See also R. J. Baxter, J. Phys. A 19 (1986) 2821. 
[21] R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 5165. 

[22] R. Shrock and S.-H. Tsai, J. Phys. A 30 (1997) 495; Phys. Rev. E55 (1997) 6791; Phys. Rev. E56 (1997) 2733. 
[23] R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 1342; Phys. Rev. E56 (1997) 3935. 
[24] R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 4111. 
[25] J. Salas and A. Sokal, J. Stat. Phys. 86 (1997) 551. 

[26] N. L. Biggs, R. M. Damerell, and D. A. Sands, J. Combin. Theory B 12 (1972) 123; N. L. Biggs and G. H. J. 
Meredith, J. Combin. Theory B 20 (1976) 5. The region diagram for the infinite-n limit of ladder graphs, shown 
as Fig. 1 of the first paper, was not completely correct; the right-most complex-conjugate curves of the boundary 
B that were thought to terminate in region Ri actually do not, but instead continue all the way in to the real 
axis, meeting sd, q — 2 and thus enclosing two further regions. See Ref. j2^ and Fig. 3 of Ref. [ pl| . 

[27] S. Beraha, J. Kahane, and N. J. Weiss, J. Comb. Theory B 28 (1980) 52. 

[28] R. C. Read and G. F. Royle, in Graph Theory, Combinatorics, and Applications (New York, Wiley, 1991), vol. 
2, p. 1009; R. C. Read, in the Proceedings of the Fifth Caribbean Conference on Combinatorics and Computing, 
Barbados, 1988. The region diagram for the infinite-n limit of ladder graphs, shown as Fig. 4 and 13 in these 
respective papers, was not completely correct; on the left, the complex-conjugate curves of the boundary B that 
were thought to terminate on the unit circle at g ~ 0.191 ±0.982i (i.e., at g = e*'*= with = arccos((3— \/5)/4)) 
actually do not, but instead continue all the way to the real axis, meeting at the origin, q = 0, and thus enclosing 
another region, as shown in Fig. 3 of Ref. [^ . Professor Read has kindly informed us (private communication) 
that he has calculated recursion relations and chromatic zeros for the Ly = 3 square strip in unpublished work. 

[29] M. Rocek, R. Shrock, and S.-H. Tsai, to appear. 

[30] J. F. Nagle, J. Combin. Theory 10 (1971) 42; G. A. Baker, Jr., J. Combin. Theory 10 (1971) 217; D. Kim and 
L G. Enting, J. Combin. Theory, B 26 (1979) 327. 



47 



[31] B. Griinbaum and G. Shephard, Tilings and Patterns (Freeman, New York, 1987). 
[32] V. Matveev and R. Shrock, J. Phys. A 28 (1995) 5235. 
[33] V. Matveev and R. Shrock, Phys. Lett. A221 (1996) 343. 

[34] S. Lcfschetz, Algebraic Geometry (Princeton Univ. Press, Princeton, 1953); R. Hartshorne, Algebraic Geometry 

(Springer, New York, 1977). 
[35] V. Matveev and R. Shroek, Phys. Lett. A204 (1995) 353. 
[36] R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 5184. 
[37] P. P. Martin, J. Phys. A 19 (1986) 3267; J. Phys. A 20 (1986) L601. 
[38] E. J. Farrell, Discrete Math. 29 (1980) 161. 



48 



